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Abstract. Real networks can be classified into two categories: fractal networks and non-fractal networks. 
Here we introduce a unifying model for the two types of networks. Our model network is governed by a 
parameter q. We obtain the topological properties of the network including the degree distribution, average 
path length, diameter, fractal dimensions, and betweenness centrality distribution, which are controlled by 
parameter q. Interestingly, we show that by adjusting q, the networks undergo a transition from fractal to 
non-fractal scalings, and exhibit a crossover from 'large' to small worlds at the same time. Our research 
may shed some light on understanding the evolution and relationships of fractal and non-fractal networks. 
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1 Introduction 

The past ten years have witnessed a considerable interest 
in characterizing and understanding the topological prop- 
erties of networked systems jTJIHEIHM] . It has been estab- 
lished that small-world property [6] and scale-free behav- 
ior [7j are the two most fundamental concepts constitut- 
ing our basic understanding of the organization of many 
natural and social systems. A serial of recent research in- 
dicate that these two features often go along, and have im- 
portant consequences on almost every aspect on dynamic 
processes taking place on networks [3,4]. The small- world 
characteristic means that the node number of network (or- 
der) increases exponentially with the average path length 
(APL), and thus leads to the general understanding that 
complex scale-free small- world networks are not invariant 
or topologically fractal, since fractal networks implies that 
there is a power-law relation between the network order 
and its APL. 

More recently, by using a renormalization procedure 
based on the box-counting method, Song, Havlin and Makse 
discovered that some real-life networks exhibit fractal scal- 
ing [8,9]. The fractal topology can be characterized via 
two relevant exponents: fractal dimension ds and degree 
exponent of the boxes dk- The fractal dimension ds is 
measured by the scaling of the minimum number of boxes 
Nb of linear size is that is needed to cover the network 
with order N, in other words, Nb/N ~ £ B ~ dE '. Similarly, 
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one can identify the degree exponent of the boxes through 
the relation kB{^B)/khub ~ ^B dfc , where fcs(€s) is the de- 
gree of each node of the renormalized network, and khub 
the maximum degree of nodes inside each box of original 
network. In fractal scale-free networks with degree distri- 
bution P{k) ~ A: -7 , the three indexes 7, ds and dk are 
related by 7 = 1 + d B /d k [HUHl • 

According to the presence of fractal scaling or not, 
networks can be assorted into two categories [5]: In the 
presence of fractal behavior, a network is said to be frac- 
tal; in contrast, if a network exhibits no fractal scaling, 
it is defined as non-fractal. Examples of the first class of 
networks include the World Wide Web (WWW), the ac- 
tor collaboration networks, metabolic networks, and yeast 
protein interaction networks. And instances of the second 
type are the Internet, and most model networks such as 
the Barabasi -Albert (BA) network [7], the Watts-Strogatz 
(WS) network [6], and the Erdos-Renyi (ER) random graph [TU] . 
In addition to different topological aspects [8lfTT] , the two 
types of networks also have distinct consequences regard- 
ing the physics of dynamical models running on them [121 

nannaiTini. 

Given the fact that real networks are either fractal or 
non-fractal, it is consequently of fundamental importance 
to understand the growth mechanisms and uncover the 
origins of different kinds of networks. To this end, a wide 
variety of models have been presented [T][21E]. However, 
previous network models, to the best of our knowledge, 
can generate either fractal networks or non-fractal ones, 
but rarely both [SJ. Thus, it seems quite natural and inter- 
esting to set up a unifying framework looking for a deeper 
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connection between fractal and non-fractal networks. This 
is the purpose of the current work. 

In this paper, by introducing a simple network growth 
process we propose a unifying scenario for fractal and non- 
fractal networks. We analytically obtain many structural 
characteristics of the network, including degree distribu- 
tion, average path length, diameter, fractal dimension, and 
betweenness centrality. The degree distribution obeys a 
power law with an exponent varying continuously. The ob- 
tained results on APL and diameter show that the network 
undergoes a transition from a small- to large-world net- 
work. More interestingly, the network exhibits a crossover 
behavior between fractality and non-fractality. 

2 Network model 

This section is devoted to network construction and com- 
putation of some related quantities. 




Fig. 1. Iterative construction method of the networks. Each 
link is replaced by either of the two paths on the right-hand side 
of the arrow with a certain probability, where each ■ stands 
for a new node. 



As will be shown in next section, when q = 1, the 
network is a small world with infinite fractal dimension, 
its average path length (APL) grows logarithmically with 
node number. In the particular case of q = 0, the network 
is a 'large' world with a finite fractal dimension of 2, its 
APL scales exponentially with network size. Except these 
two limiting cases of q = 1 and q = 0, for other q (0 < 
q < 1), the networks are growing stochastically. Varying 
q in the interval (0,1), the networks exhibit a transition 
from small to large world, and simultaneously show an 
interesting phenomenon with a transition from non-fractal 
to fractal behavior. 

Note that the current model is similar to the prob- 
abilistic model introduced by Song, Havlin, and Makse 
(SHM Model) in [5]. The construction step with probabil- 
ity q in the current manuscript is analogous to the Mode I 
growth process in the SHM Model, involving strong hub- 
hub attraction. And the step with probability 1 — q is anal- 
ogous to the Mode II growth process in the SHM Model, 
involving hub-hub repulsion. The probabilistic mixing of 
these two types of construction steps leads to the same 
results in both of the SHM Model and the current model: 
when g=l (or pure Mode I) is involved, the network 
has infinite fractal dimension and shows small- world scal- 
ing; any q < 1 (or mixing of Modes I and II) leads to 
a finite fractal dimension. However, in [9] only part of 
the properties was addressed. Here we will present an ex- 
haustive analysis of various characteristics, including addi- 
tional features (such as average path length and between- 
ness centrality) not calculated in [9|. 



2.2 Order and size 



2.1 Construction algorithm 

The proposed evolving networks (graphs) have a treelike 
structure, which are constructed in an iterative way as 
shown in Fig. Q] Let T t {t > 0) denote the networks after 
t iterations. For t = 0, the networks growth begins from 
two nodes (vertices) connected by an edge (link or bond) . 
For t > 1, Tt is obtained from T t -\. We replace each edge 
in T t _i either by the path on the top right of Fig. [1] with 
probability q, or by the path on bottom right with comple- 
mentary probability 1 — q. In other words, T t is generated 
from T t _i by adding k v leaves to each node v, where k v is 
the degree of node v in T t _i. Then, the resulting graph is 
further modified by expanding each edge with probability 
1 — q, which was already present in T t -\. Expansion of 
an edge (uv) means: removal of (uv) and creation of an 
additional node w with edges (mw) and (vw). (Note that 
real systems may exhibit such an evolving mechanism. For 
example, this mechanism has been used to model the net- 
work evolution of connected minima on a potential energy 
landscape [15].) This procedure is iterated t times, with 
the needed graphs obtained in the limit t — > oo. 



Now we compute some related quantities such as the num- 
ber of total nodes and edges in T t , called network order 
and size, respectively. It should be mentioned that as q is a 
real number, we assume that all variables concerned with 
q change continuously. Note that a similar assumption was 
used in Refs. [1,2,3,4], which is valid for large t. Let L v (t) 
be the number of nodes generated at step t, E t the total 
number of edges present at step t. Then L v (0) = 2 and 
.Eq = 1. By construction (see Fig. [T|), at each time step, 
each existing edge is replaced either by three edges with 
probability q or by four edges with complementary proba- 
bility 1-q. Thus, S t = qx3E t -i + (l-q)x4E t _i = (4-g)' 
(t > 0). At the same time, each existing edge yields two 
or three new nodes with probability q or 1 — q, this leads 
to L v {t) =qx 2E t -! + (1 - q) x 3E t -i = (3 - ?)(4 - qf" 1 
(t > 1). Then the number of total nodes N t present at 
step t is 

t 

N t = Y J L v{U) = {±-q) t + l. (1) 

The average node degree after t iterations is (k) t = 
which approaches 2 for large t, coinciding with the treelike 
structure of networks. 
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3 Topological properties 

In this section, we will show that the tunable parameter q 
in the above construction algorithm controls the relevant 
features of the networks. 



3.1 Degree distribution 

When a new node i is added to the networks at a certain 
step t{ (ti > 1), it has a degree of either 1 or 2. We denote 
by ki(t) the degree of node i at time t. By construction, 
the degree fcj(i) evolves with time as fc,(t) = 2 ki(t — 1). 
That is to say, the degree of node i increases by a factor 
2 at each time step. Thus, the degree spectrum of the 
networks is discrete. In network T t all possible degree of 
nodes is 1, 2, 2 2 2 3 , . . ., 2*" 1 , 2*. 

Let L\(ti) and L^{ti) be the separate number of new 
nodes with degree 1 and 2 that were born at step U. Ac- 
cording to the construction algorithm, we have Life) = 
q x 2E ti -i + (1 - g) X 2E ti -i = 2 (4- q) 1 -" 1 and L 2 (U) = 
(1 - q)E U -i = (1 - g)(4 - g)* 4-1 . Then in network T t , 
the expected number of nodes of degree k — 2 l ~ m is 
n k = L x {m) + L 2 (m + 1) = (6 - 5g+ g 2 )(4 - g)™- 1 . 

Since the degree spectrum of the networks is not con- 
tinuous. It follows that the cumulative degree distribu- 
tion [3] is given by P cum (k) = where A tjfe = X) fe ,> fe n k> 
is the number of nodes whose degree is not less than k. 
When t is large enough, we find P cum (fc) w ^.-ln(4- 9 )/in2_ 
So the degree distribution P(k) of the networks follows a 
power-law form P(k) *~ fc -7 with the exponent 



which is approximately equal to 



13x2* 



7 



1 



ln(4 - q) 
lri2 : 



(2) 



which is a monotonically decreasing function of q. As g 
increases from to 1, 7 drops form 2 to 1 + |"-| [TT?] . 



3.2 Average path length 

Shortest paths play an important role both in the trans- 
port and communication within a network and in the char- 
acterization of the internal structure of the network [5]. 
Let dij represent the shortest path length from node i to 
j, then the average path length (APL) d t of T t is defined 
as the mean of dij over all couples of nodes in the network, 
and the maximum value D t of dij is called the diameter 
of the network. APL is relevant in many fields regarding 
real-life networks and has received much attention 20J . 

For general q, it is difficult to derive a closed formula 
for the APL dt of network T t . But for two limiting cases of 
q — and q = 1, both the networks are deterministic ones, 
which allows one to obtain the analytic solutions for APL. 
The detailed exact derivation about APL is included in 
the Appendix section. The obtained results show that the 
APL for these two particular cases presents qualitatively 
disparate behaviors: For q = 0, it is found that 



+ 7 x 4* + 13 x 
14(1 + 4*) 



(3) 



Ni 



4* for large network, so dt 



14 

A, 



1/2 



for large t. Since 
indicating that 



dt grows as a square power of the network order N t . This 
phenomenon is similar to that of the two-dimensional reg- 
ular lattice [21] . Thus, the network corresponding to q = 
is not a small world. For q = 1, we find 



d t = 



2(1 + 2 x 3* +t x 3*) 
3(1 + 3') ' 



(4) 



which approximates m the infinite t, implying that 
the APL shows a logarithmic scaling with network order. 
Therefore, in the specific case of q = 1, the network ex- 
hibits a small- world behavior. 

Thus, when we tune q from to 1, the networks un- 
dergo a transition from a 'large' to small world. We stress 
that such a transition has already been observed in some 
previously studied models |22j . 



3.3 Diameter 

As mentioned in preceding subsection, the diameter of a 
network is defined as the longest shortest path between 
all pairs of nodes, characterizing the maximum commu- 
nication delay in the network. Although we do not give 
a closed formula of APL of T t for general q, here we will 
provide the exact result of the diameter of T t denoted by 
D t . 

We first address the network of q = 1 case, where the 
shortest distances between existing node pairs are not al- 
tered when new nodes enter the systems. For this par- 
ticular case, Do = 1. At each time step, the diameter 
of the network increases by 2. Then the diameter of T t 
is D t — 2 1 + 1 and thus scales logarithmically with the 
network order, showing a similar behavior as that of the 
average path length. Since small diameter is consistent 
with the concept of 'small- world', the additive growth in 
the diameter with time also (as the APL) suggests that 
the network for q = 1 case is a small world. 

For < q < 1, the addition of new nodes affects fun- 
damentally the distances between existing node pairs. By 
construction algorithm, for any existing couple of nodes 
connected by a link, after a generation of evolution, the 
distant I between this pair of nodes may be equal to 1 or 
2 with probability q and 1 — q, respectively. Thus the ex- 
pected value of I is 2 — q. Using this result and considering 
that the networks are treelike, we can derive the following 
recursive relation for expected D t : 



D t = (2-g)A-i+2. 



(5) 



Since Do = 1 , we can resolve Eq. ([5]) to obtain the average 
of network diameter as 



Dt 



l-g' 



(G) 



which grows as a power of time t. 



4 



Zhongzhi Zhang et al.: Transition from fractal to non- fractal scalings in growing scale- free networks 



3.4 Fractal dimension 

To determine the fractal dimension, we distinguish two 
cases: < q < 1 and q = 1. In the case of < q < 
1, we follow the mathematical framework introduced in 
Ref. [9j. We are concerned about three quantities: network 
order N t , network diameter D t , and degree k u (t) of a given 
node u. By construction, we can easily see that in the 
infinite t limit, these quantities grow obeying the following 
relations: N t ~ (4 - q) N t - U D t ~ (2 - q) D t -x, k u (t) = 
2k u (t — 1). Thus, for large networks, Nt, k u (t) and D t 
increase by a factor of Jn = 4 — q, fk — 2, and fr> = 2 — q, 
respectively. 

From above obtained microscopic parameters demon- 
strating the mechanism for network growth, we can derive 
the scaling exponents: the fractal dimension cLb = = 

}"|^~^ and the degree exponent of boxes dk — j^y^ — 

1,1 2 According to the scaling relation of fractal scale- 



free networks, the exponent of the degree distribution sat- 
isfies 7 = l + ^ a - = l + ' §i vm g t ne same 7 as that 
obtained in the direct calculation of the degree distribu- 
tion, see Eq. @. 

For q = 1, although the number of their nodes in- 
creases exponentially, its diameter grows linearly with time. 
Thus, in this case the network has infinite dimension and 
does not present a fractal topology. 



Thus, the BC of these d(g) nodes is not less than g. On 
the other hand, in the whole network there are N t /g such 
areas, each of which includes g nodes. Then, the total 
number of nodes with BC b not less than g in the whole 
network is 

nib > g) ~ d(g) x — ~ N t x (8) 
9 

Therefore, the cumulative BC distribution is 



Pcum (^) 



n(b > g) 
N t 



-(l-i/ds) 



which implies P(b) ~ b 76 with 

76 = 2--L 

ds 



(9) 



(10) 



Thus, 7b is a increasing function of fractal dimension ds, 
the larger the fractal dimension ds, the larger the ex- 
ponent 7b. Eq. (fTOf shows that 7b = 2 holds only for the 
non- fractal networks with ds — * 00 [28] , where a relatively 
small number of hub nodes bear large BC. In contrast, in 
fractal networks, the BC of a lot of 'small' nodes can be 
compared with that of hub nodes. 



4 Conclusions 



3.5 Betweenness centrality 

Bctweenness centrality (BC) of a given node is the ac- 
cumulated fraction of the total number of shortest paths 
going through the node over all node pairs [24jl23], More 
precisely, the betweenness of a node i is 



(7jk 



(7) 



where Ujk is the total number of shortest path between 
node j and fc, and <Jjk{i) is the number of shortest path 
running through node i. 

We now investigate the BC distribution of nodes. It 
was claimed in previous studies that for scale-free net- 
works, the BC distribution P(b) of nodes obeys a power 
law with the exponent 7b = 2, i.e., P(b) ~ b~ 2 . It was 
also suggested that the exponent 7b = 2 is universal for 
all scale-free networks [25]. Next, we will present that the 
exponent 7b is not universal and constant (at least for 
fractal trees), but varies significantly as a function of ds 
(or a function or q). 

In order to obtain the exponent 7b of BC distribution 
of our networks, we resort to a heuristic argument similar 
to that applied in [55] . Notice that all of our networks are 
treelike and fractal (except the case of q = 1). Consider- 
ing one network with dimension ds, for a small area of 
the network containing g nodes, its average path length is 
typically d(g) — g 1 ^ 3 [27]. All nodes in this small region 
can reach the rest nodes in the network via d(g) nodes. 



In the paper, by introducing a parameter q, we have pre- 
sented a simple network growth process to generate a uni- 
fied model for fractal and non-fractal networks. This pro- 
cess was shown to lead to a rich behavior for the network 
structure. Various relevant topological properties have been 
determined depending on the model parameter q. It has 
been shown that both degree and betweenness centrality 
distributions of nodes have a power-law tail and the char- 
acteristic exponents change continuously with the param- 
eter q. Other structural properties including the APL, di- 
ameter, and fractal dimension have been obtained as well, 
which indicate that the model undergoes a crossover from 
large-world to small-world networks, and simultaneously 
exhibits a transition from fractal to non-fractal behaviors. 

In spite of the simplicity, our minimal model can cap- 
ture the essential characteristics and correlations of fractal 
and non-fractal networks. It is helpful for understanding 
the growth mechanisms and evolutions of the two different 
sorts of networks. Finally, it should be mentioned that, al- 
though we only studied treelike networks, in a similar way, 
one can construct models considering the effect of loops, 
whose general properties are similar to those of the model 
investigated in the present work. 
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Appendix A: Derivation of the average path 
length for two limited cases 

Following an algebraic method similar to that introduced 
in p~2] , we can compute the average path length ( APL) for 
two limiting deterministic cases. By definition, the APL 
for T t is defined as 



d t = 



N t (N t - l)/2 ' 



where 



S t = 



E 



(A.1) 



(A.2) 



denotes the sum of the shortest path length between two 
nodes over all pairs. For the two particular cases q = 
and ?=1, Both of the networks have a self-similar struc- 
ture that allows one to calculate dt analytically. The self- 
similar structure is obvious from an equivalent network 
construction method: to obtain Tt+i, one can make some 
copies of T t and join them in the hub nodes. 



1 The q = case 

As shown in Fig. for the q = case, the network T t +i 
may be obtained by the juxtaposition of four copies of T t , 
which are labeled as T? , 9 = 1,2,3,4. Then we can write 



the sum St+i as 



4 S t + A t 



(A.3) 



where At is the sum over all shortest paths whose end- 
points are not in the same F t branch. The solution of 
Eq. (|A~3|) is 



A 1 - 1 Si 



t-i 



~ l A x . 



(A.4) 



The paths that contribute to A t must all go through at 
least one of the three edge nodes (i.e., E, F and G in 
Fig. [Ub)) at which the different T t branches are con- 
nected. The analytical expression for A t , called the cross- 
ing paths, is found below. 

Denote A"' 13 as the sum of all shortest paths with end- 
points in T" and Tf . If T" and Tf meet at an edge node, 
A^ rules out the paths where either endpoint is that 



shared edge node. If T" and Tf do not meet, Af'* 3 ex- 
cludes the paths where either endpoint is any edge node. 
Then the total sum A t is 



At = A]' 2 + A 1 / + Al' 4 + A 2 t z ' * 2A ' ^ 



A^ + AY. (A.5) 



By symmetry, A)? = Af A = A^ A and A\' A = Z\f 4 , 
so that 



A t = 3A 



1.2 



2Al> 3 + A t 



1,4 



In order to find A t ' , A t ' , and A t ' , we define 

st= ^ d i,E ■ 
ieT^.i^E 



(A.6) 



(A.7) 



Considering the self-similar network structure, we can eas- 
ily know that at time t + 1, the quantity st+i evolves re- 
cursively as 

s t +i = 2s t + [s t + 2\N t - 1)] + [s t + 2 t+1 {N t - 1)] 
= 4s t + 3x2 3 *. (A.8) 

Using si = 7, we have 

S( = 2 2 " + 3x2" (A.9) 

On the other hand, by definition given above, we have 



AY = 



E 



(diE + dj E ) 



iJ¥=E 



{N t - 1) E d * E + ( N * ~ J ) E d J E 



i=£E 

2(N t - 1) diE 

iET} i^E 

2{N t -l)s t , 



(A.10) 



A 



1.3 



and 



/ - E di J 

i£T?,i^E 

= ^ (diE + d E F + d jF ) 

i£Tt,i^E 

= 2(iV t -l) S( + (JV t -l) 2 x2 t , (A.ll) 



A t A = E di i 

ieTt,i^E 

= E (diE + d E G + d jG ) 

i£Tt,i^E 

= 2(JV t -l) Si + (JV ( -l) 2 x2 i+1 , (A.12) 
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Fig. 2. (Color online) The network growth process for the particular case q — 0. (a) Illustration of the first four evolution steps, 
(b) Second construction method of the network for q — case that highlights self-similarity: The graph after t + 1 construction 
steps, T t +i, is composed of four copies of Tt denoted as Tf (8 = 1, 2, 3, 4), which are connected to one another as above. 



where <1ef = 2 l and cLeg = 2 t+1 have been used. Substi- 
tuting Eqs. COB, fOTj) and (TA~T2j) into Eq. (fO|l . we 
obtain 

A - 12(JV t - 1) s t + 4 (N t - If x 2* 

= 13 x 32* + 3 x 16*. (A.13) 

Inserting Eqs. (|A. 13|) for A x into Eq. (|A.4|1 . and using 
5i = 20, we have 

4* 

S t = — x (8 + 7 x 4* + 13 x 8*) . (A.14) 
28 

Inserting Eq. (|A.14|) into Eq. (|A. 1|) . one can obtain the 
analytical expression for dt in Eq. ([3]). 

2 The q = 1 case 

Using analogous analysis for the q — case, we can calcu- 
late the APL dt for the case of q = 1. For simplicity, we use 
the same symbols used in last subsection to represent the 
identical notions. For q — 1 case, the network Tt+i may 
be obtained by joining at the hubs (the most connected 

nodes) three copies of T t labeled as ip = 1,2,3 

see Fig. [3] Then one can write the sum over all shortest 

paths St+i as 

S t+1 =3S t + A t . (A.15) 
The solution of Eq. (|A.15P is 

t-i 

S t = 3*- 1 1 S , i + ^3*- T -M r . (A.16) 

r=l 

The paths contributing to Z\t must all go through at 
least either of the two hubs (X and Y) where the three 
different T t branches are joined. The crossing paths A t is 
given by 

A t = A]' 2 + A 2 ' 3 + A]' 3 , (A.17) 



where A\ 2 rules out the paths where either endpoint is 
node X, A t ' rules out the paths where either endpoint is 
node Y, and A t ' excludes the paths with an endpoint is 
either X or Y. Again by symmetry, A t ' — A t ' , so that 

A t = 2A]' 2 + A]' 3 . (A. 18) 

In this case, the quantity St+i evolves as 

s t+ i = 2 at + [at + (N t - 1)] 

= 3s t + 3 t . (A.19) 

Since a% = 4, Eq. (|A.19j) is solved inductively: 

s t = (t + 3) X 3*" 1 . (A.20) 

similarly, 

A 1 ' 2 = E (dix + djx) 

= 2(Nt-l)s u (A.21) 

and 

A]' 3 = E ( d iX + dxv + d jY ) 

= 2{N t - 1) s t + (N t - l) 2 . (A.22) 

Substituting the obtained expressions in Eqs. (|A.21[) and 
(|A.22|) into Eq. (|A.18[) . the crossing paths A t is found to 
be 

A t = 7 x 9* + 2t x 9*. (A.23) 

Inserting Eq. (|A.23[) into Eq. (|A.16j) and using the initial 
condition Si = 10, we have 

S t = 3- 1+t (1 + 2 x 3* + t x 3*) . (A.24) 

Substituting Eq. (|A.24|) into (|A.1|) . the exact expression 
for the average path length is obtained as shown in Eq. (0J . 



Zhongzhi Zhang et al.: Transition from fractal to non- fractal scalings in growing scale- free networks 



7 



• • ®— • — ®- 

t=0 t=l 




t=2 

t=3 



(b) 




Fig. 3. (Color online) The network growth process for the specific case q = 1. (a) Scheme of the first four iterative processes, 
(b) Second construction means of the network for q = 1 case: The graph after t + 1 construction steps, Tt+i, can be obtained 
by joining three copies of Tt denoted as Tf (ip — 1, 2, 3), which are connected to each other at the two edge nodes (X and Y). 



References 

1. R. Albert and A.-L. Barabasi, Rev. Mod. Phys. 74, 47 
(2002). 

2. S.N. Dorogovtsev and J.F.F. Mendes, Adv. Phys. 51, 1079 
(2002). 

3. M.E.J. Newman, SIAM Rev. 45, 167 (2003). 

4. S. Boccaletti, V. Latora, Y. Moreno, M. Chavez and D.-U. 
Hwanga, Phys. Rep. 424, 175 (2006). 

5. L. da. F. Costa, FA. Rodrigues, G. Travieso, and P.R.V. 
Boas, Adv. Phys. 56, 167 (2007). 

6. D.J. Watts and H. Strogatz, Nature (London) 393, 440 
(1998). 

7. A.-L. Barabasi and R. Albert, Science 286, 509 (1999). 

8. C. Song, S. Havlin, H. A. Makse, Nature 433, 392 (2005). 

9. C. Song, S. Havlin, H. A. Makse, Nature Phys. 2, 275 
(2006). 

10. P. Erdos, A. Renyi, Pub. Math. Insti. Hung. Acad. Sci. 5 
(1960) 17. 

11. M. Kitsak, S. Havlin, G. Paul, M. Riccaboni, F. Pammolli, 
and H. E. Stanley, Phys. Rev. E 75, 056115 (2007). 

12. M. Hinczewski and A. N. Berker, Phys. Rev. E 73, 066126 
(2006). 

13. Z. Z. Zhang, S. G. Zhou, and T. Zou, Eur. Phys. J. B 56, 
259 (2007). 

14. Z. Z. Zhang, S. G. Zhou, T. Zou, and J. H. Guan, (unpub- 
lished) . 

15. H. D. Rozenfeld, S. Havlin, and D. ben-Avraham, New J. 
Phys. 9, 175 (2007). 

16. H. D. Rozenfeld and D. ben-Avraham, Phys. Rev. E 75, 
061102 (2007). 

17. M. Hinczewski, Phys. Rev. E 75, 061104 (2007). 

18. C. P. Massen and J. P. K. Doye, J. Chem. Phys. 127, 
114306 (2007). 

19. S. Jung, S. Kim, and B. Kahng, Phys. Rev. E 65, 056101 
(2002). 

20. A. Fronczak, P. Fronczak, and J. A. Holyst, Phys. Rev. E 
70, 056110 (2004); J. A. Holyst, J. Sienkiewicz, A. Fron- 
czak, P. Fronczak, and K. Suchecki, Phys. Rev. E 72, 



026108 (2005); S. N. Dorogovtsev, J. F. F. Mendes, and J. 

G. Oliveira, Phys. Rev. E 73, 056122 (2006); Z.Z. Zhang, 
L.C. Chen, S.G. Zhou, L.J. Fang, J.H. Guan, and T. Zou, 
Phys. Rev. E 77, 017102 (2008). 

21. M. E. J. Newman, J. Stat. Phys. 101 819 (2000). 

22. K. Klemm and V. M. Egufiuz, Phys. Rev. E 65, 057102 
(2002); L. Tian, C.-P. Zhu, D.-N. Shi, Z.-M. Gu, and T. 
Zhou, Phys. Rev. E 74, 046103 (2006); Z. Z. Zhang, S. G. 
Zhou, Z. Shen, and J. H. Guan, Physica A 385, 765 (2007); 
Z. Z. Zhang, S. G. Zhou, Z. Y. Wang, and Z. Shen, J. Phys. 
A 40, 18863 (2007); S. N. Dorogovtsev, P. L. Krapivsky, 
and J. F. F. Mendes, Europhys. Lett. 81, 30004 (2008). 

23. M. E. J. Newman, Phys. Rev. E 64, 016132 (2001). 

24. C. L. Freeman, Sociometry 40, 35 (1977). 

25. K.-I. Goh, B. Kahng, and D. Kim, Phys. Rev. Lett. 87, 
278701 (2001). 

26. L. A. Braunstein, Z. Wu, Y. Chen, S. V. Buldyrev, T. 
Kalisky, S. Sreenivasan, R. Cohen, E. Lopez, S. Havlin, and 

H. E. Stanley, Int. J. Bifurcation Chaos 17, 2215 (2007). 

27. Fractals in Science, edited by A. Bunde and S. Havlin 
(Springer, Berlin, 1996). 

28. C.-M. Ghima, E. Oh, K.-I. Goh, B. Kahng, and D. Kim, 
Eur. Phys. J. B 38, 193 (2004); Z. Z. Zhang, S. G. Zhou, 
L. C. Chen, J. H. Guan, L. J. Fang, and Y. C. Zhang, Eur. 
Phys. J. B 59, 99 (2007). 

29. E. Bollt, D. ben-Avraham, New J. Phys. 7, 26 (2005); F. 
Cornelias, H. D. Rozenfeld, D. ben-Avraham, Phys. Rev. 
E 72, 046142 (2005). 



